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We show that the dark energy effects can be modeled by using an Ismg perfect fluid with net- 
work interactions, whose low redshift equation of state, i.e. uq, becomes ljo — —1 as in the ACDM 
model. In our picture, dark energy is characterized by a barotropic fluid on a lattice in the equi- 
librium configuration. Thus, mimicking the spin interaction by replacing the spin variable with 
an occupational number, the pressure naturally becomes negative. We find that the corresponding 
equation of state mimics the effects of a variable dark energy term, whose limiting case reduces to 
the cosmological constant A. This permits us to avoid the introduction of a vacuum energy as dark 
energy source by hand, alleviating the coincidence and fine tuning problems. We find fairly good 
cosmological constraints, by performing three tests with supernovae la, baryonic acoustic oscillation 
and cosmic microwave background measurements. Finally, we perform the AIC and BIG selection 
criteria, showing that our model is statistically favored with respect to the GPL parametrization. 
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I. INTRODUCTION 

Recent cosmological observations pointed out that the 
universe is undergoing an accelerated expansion [l|, y]- 
Unfortunately, the physical mechanism which drives the 
observed cosmic speed up is so far unclear. Moreover, 
standard pressureless baryonic matter is inadequate by 
itself to characterize the universe acceleration, even by 
assuming the additional presence of cold dark matter [3|- 
Igj. As a standard landscape, it is possible to postulate 
the existence of a further ingredient, dubbed dark en- 
ergy (DE) [7h11|. Dark energy behaves as a weakly in- 
teracting anti-gravitational fluid, described by a negative 
equation of state (EoS). Even though the nature of such 
a fluid has not been clarified, a wide number of differ- 
ent paradigms have followed each other, spanning from 
slowly rolling scalar field, known as quintessence [l2l [l3| , 
lattice of topological defeats [IJ, to barotropic fluids 
|15| or modifications of Einstein gravity |16l - KL8| and so 
forth (for further details see [13 and references therein). 
One of the simplest way to explain the observed cosmic 
speed up is provided by the so-called ACDM paradigm. 
In particular, the ACDM model leads to the introduc- 
tion of a vacuum energy cosmological constant A, and 
assumes a total matter content, rim, given by the sum of 
baryonic and dark matter densities [20, [2l| . The corre- 
sponding EoS is constant as the universe expands, lead- 
ing moreover to a negative and constant pressure j22| . 
Even though the model is in a fairly good agreement 
with current observations, it suffers from two profound 
shortcomings, i.e. the fine tuning and coincidence prob- 
lems [13, 19, 2j,[2J]. To alleviate these two issues, the 
DE density, p, and its corresponding pressure, 'P{p), are 
thought to evolve separately with significative departures 



from standard matter. Thus, possible extensions of the 
ACDM model are frequently characterized by assuming a 
time-variable EoS, i-o{z) = — ^, evolving as — 1 < a; < 0, 
at z = 0. 

In this work, we investigate how to model DE by con- 
side ring an Ising network-interacting fluid on a lattice 
[13 . |25| - [27| . We show that a fiuid, whose interaction is 
provided by a series of networks, rnay predict effects due 
to a negative pressure, at z <C 1 [23 ■ Thence, the DE na- 
ture may be approximated by a network interacting Ising 
fiuid, simply postulating the validity of standard thermo- 
dynamics, in the equilibrium configuration. In particu- 
lar, one can assume that DE is modeled in analogy to 
the case of Ising chains. In our model, the role played 
by spin variables is replaced by occupational numbers. 
In doing so, one mimics the spin interaction without the 
need of real spin chains. The main advantage is to re- 
cover the physical properties of a spin interacting fluid 
in cosmology, predicting repulsive effects, able to explain 
the observed late time DE. The procedure of mimicking 
an Ising system has been extensively investigated in the 
literature. Examples have been considered in the case of 
liquid-gas phase transitions for atomic systems [29[. 

This paper is structured as follows. In Sec. II, we 
show the main features of our model. We focus on its 
thermodynamic interpretation, with particular attention 
to relate these results with the Friedmann equations. In 
Sec. Ill, we derive the cosmological model and the cor- 
responding EoS. We investigate the acceleration param- 
eter and its variation. We find out theoretical limits on 
the observable quantities and we determine the transition 
redshift, at which the acceleration starts. In Sec. IV, 
we perform three experimental tests, with supernovae la 
(SNela), baryonic acoustic oscillation (BAO) and cosmic 



microwave background (CMB), in order to constrain our 
model. In Sec. V, we adopt the AIC and BIG selection 
criteria, for comparing our model with alternative ap- 
proaches. Finally, Sec. VI is devoted to conclusion and 
perspectives of our work. 



II. THE ISING FLUID WITH NETWORK 
INTERACTIONS 

In this section, we describe the cosmological conse- 
quences of assuming an Ising fluid with network inter- 
actions, as a source of DE. Let us first notice that in 
standard lattice models, the particle description is for- 
mally represented by assuming a grid with an occupa- 
tional variable Ci. Its value is zero, if the site is empty, 
and one if occupied by particles. This simple picture is 
also known in the literature as "bit" gas model, in anal- 
ogy to computational science [30, l3l| . Its use is usually 
adopted in different fields of physics, especially in order to 
model extended sites by numerical computations through 
Monte Carlo simulations, ranging from condensed matter 
to quantum computing and particle physics |32l440l | . Un- 
der these hypotheses, the occupational variables Ui are 
actually analogous to spin variables in the well known 
Ising model [4ll |42| . However, the physical meaning of 
Ui is basically different from Ising spins. In particular, 
Ui is not an intrinsic physical property of the system, 
but only a way to discriminate the presence and absence 
of particles. Moreover, once the particle interaction is 
specified, our lattice fluid may lead to a negative EoS, 
in particular regions of the phase space [4j, |4J] . We will 
consider this property, in order to describe the DE effects 
in a homogeneous and isotropic universe. 



and temperature are functions of the redshift, we choose 
an unitary lattice spacing and a hard sphere interaction, 
which consists in introducing an excluded volume by pro- 
hibiting the multiple occupancy of particles on a given 
lattice site [46|. In order to obtain the EoS of our model, 
we can take into account a finite region of volume V at a 
given temperature T, where TV particles are confined to 
move only on d-dimensional discrete lattice points, inside 
the region under interest [43]. The volume is expressed 
by V = ft. • L''~^ with h the height and L'^~^ the area of 
the side wall. Thus, the entropy of the system becomes 



S = kB\nQ.{M,V). 



(3) 



where Q.{J\f ,V) represents the total permutations of 
putting M particles in a volume V. In the absence of 
excluded volume interactions, the entropy of the system 
is factorized 



S = kB In n{N, V) = kB ln(a;f ) = kBU\n{uj{) , (4) 



where £7 



,A^ 



, and LOi represents the single wall proba- 



bility associated to the i—ih particle in a unitary volume. 
In the simplest case of one single particle, the entropy re- 
duces to fcslnw. In classical statistical mechanics, the 

partition function is in general Z = exp ( — — S^'H] and 

reduces to a simple form when the Hamiltonian is not 
a functional. Note that, the contributions to the parti- 
tion function come from two sources: one is associated 
to the kinetic energy of the system, while the other one 
is represented by the total potential energy. 



B. The equilibrium configuration 



A. The entropy representation 

In the Hamiltonian formalism, we can write down the 
Hamiltonian of our model, by defining a chemical poten- 
tial /i. We have 



H 



} JijCTiCTj — 2_^ l^^i 



(1) 



^,0 



where /z > in the case of micro-canonical ensem- 
ble. Here, the interaction strength between occupational 
numbers Ui is given by J^ j > 0. The Hamiltonian of Eq. 
([1]) is formally analogous to the Ising Hamiltonian in the 
mean field approximation, with the substitution [45| 



S, + l 



(2) 



where Si = ±1 represents the Ising spin variable. Once 
the Hamiltonian is known, it is easy to get the parti- 
tion function Z and to evaluate the thermodynamical 
variables. The information of network interactions is 
contained in Jij. Hence, by assuming that the volume 



In the equilibrium configuration, the contribution com- 
ing from the kinetic energy is decoupled from the con- 
figurational statistics |48l |. Hence, the equilibrium pres- 
sure can be determined by considering only the configu- 
rational properties of the system. In the case of a lattice 
configuration, the free energy is computed exactly, be- 
cause the hard sphere interaction is formally equivalent 
to introducing an excluded volume [49|. Hence, by pro- 
hibiting multiple particle occupancy on a given lattice 
site [42|, we get 



n{M, V) 



V! 



m{v-M)\ 



(5) 



By assuming that the DE density is p — ^Pa, making 
use of the Stirling approximation, we infer 



>S«-fcBV[plnp+(l-p)ln(l-p)] 



(6) 



where we plugged into Eq. @, the definition of n(7\A, V), 
i.e. Eq. ([5]). Moreover, the scaling ruler pA, is associated 
to the minimal size of lattice sites. For the sake of clear- 
ness, Pa determines whether the approximation of lattice 
fiuid holds. Thus, the term pA defines a border constant 



of the equilibrium configuration, whose physical mean- 
ing has nothing to do with the number density of lattice 
sites. For our purposes, it is convenient to assume here- 
after p\ = 1. We replace p\, in the incoming sections. 
Moreover, the pressure of our fluid, i.e. V = ujp, is given 
by 



(99 
V^-T^^-kBTHl-p)., 



(7) 



where w represents the EoS of our model. As shown in 
[501, the same results could be found in the context of 
the grand partition ensemble, without losing generality. 
The expressions of the pressure and density are in- 
terpreted in our picture as sources of DE. They enter 
the energy momentum tensor, in addition to standard 
pressureless matter. Thus, for a spatially flat homoge- 
neous and isotropic universe, i.e. ds^ — dt^ — a{t)^{dr'^ + 
j- r^d9^), the Friedmann equations are 



9 ■ 2 

r sm 



H' = 



SttG 



H + H^ 



3 

47rG 



Pt, 



(8) 



(SV + Pt) 



where both the total pressure and density, respectively 
V and Pt, are composed by pressureless matter, i.e. p„i, 
and by DE densitjo, i.e. p. 



III. COSMOLOGICAL CONSEQUENCES OF 
NETWORK INTERACTING ISING FLUID 

In this section, we are interested in finding the EoS of 
DE, modeling it as a network interacting fluid through 
the use of standard thermodynamic recipe. Hence, per- 
fect fluids are modeled by an energy momentum tensor 
of the form 



-ial3 



r"p = {p + r)u°'uP - vg 



«„/9 



,a/3 



(9) 



from that it naturally follows the conservation laws for 
energy and particle number densities. In particular, by 
assuming the redshift definition in terms of the cosmic 
time and H[z), i.e. 



§.-,i + .w=,^ 



(10) 



we get, from Eqs. ([9|) and (jlOp . the continuity equation 
for p in terms of z 



dp ^ J> + p 

dz 1 + z 



(11) 



^ Notice that the total pressure V is the pressure of DE, because 
standard matter is supposed to be pressureless. 



From Eq. ((TT|) . one gets the functional form of p 

l+uj{z) 



p oc exp 



l + z 



-dz 



and, in addition, it is possible to show that 

uj{z) 



T oc exp 



l + z 



-dz 



(12) 



(13) 



which represents the DE temperature. It is remarkable 
to notice that possible temperature measurements at dif- 
ferent epochs of the universe could discriminate whether 
DE dominates over other species or not. Current ob- 
servations seem to indicate negligible departures from 
present temperature |5l| . This does not permit us to 
fix constraints on the DE temperature today, and on 
its evolution. According to Eqs. (|12p and P^ . in the 
small redshift approximation, one recovers the equiparti- 
tion principle between p and T, having T ~ ap, with a 
a constant to be determined. Plugging the temperature 
definition in terms of p into the Eq. ([7]) , we get 



V = ap\og[l - p] . 



(14) 



Once the degrees of freedom are fixed, the correspond- 
ing cosmological model depends on the today DE value. 
From Eq. (jlip , we obtain the continuity equation for DE 



(1 + z) dp 
3 dz 



— p = ap log 



1-^ 

PA 



(15) 



Dynamics of DE Source of DE 

where we restored the definition of p\. In Eq. (|15|) . we 
split the dynamics of DE in terms of the redshift z (left 
side), from the DE source (right side), as determined by 
our model. We evaluate the Hubble rate, by expanding 
around p\ and we get 



H{z)^HoJn„^{l + z)^ + 



a 
PA 



(l + z)= 



with T an integration constant, given by 

^ ct 1 
T = -— . 

PA 1 - ^m 



(16) 



(17) 



It is worth noticing that the error due to the approxi- 
mations made in Eqs. (|6]) and (fT6)) is negligibly small. 
The numerical solution, corresponding to our logarithmic 
correction of the pressure, well approximates the exact 
solution, with an error < 10%, at small z. Thus, the EoS 
of DE can be rewritten as 



w(z) = 

whose value today is 

Wo = 



PA 



PA 



(T+ip 



PA 



(l-^n) 



(18) 



(19) 



which becomes ujq ^ —1 when 



From the 



PA ^ i-n„_ 

former relation for wq, it is evident that our model re- 
duces to ACDM as a limiting case, showing that the role 
played by the cosmological constant is actually mimicked 
through the term oc -^. To guarantee that the universe 

° PA *= 

is accelerating today, we should constrain the ratio -^. 
The degeneracy between a and pA is alleviated by the 
fact that one needs to evaluate the ratio -^, instead of 

PA ' 

a and p\ separately. Using the Friedmann equations, 
i.e. Eqs. (jH]), by keeping in mind the definition of the 
acceleration parameter: 



we find for our model 



(1 + z) dH 
H 'dz' 



q = 



3arn{l + zf-3T^uj{z 



2f]™(l + z)6 + 2 



PAJl + zY 



oj{z) 



- 1 



(20) 



(21) 



The cosmological acceleration starts at the transition red- 
shift Zacc, i-C. 



PA + 3{l-n^)^a 



9{l-n^)^a l-(l-f]^)^ 



(22) 



which approximatively reads Zacc — 0.745, when a/ p\ « 
1.341, with the indicative value ftm = 0.3. Equation ((22|) 
corresponds to the case q = 0, and leads to the epoch 
in which DE dominates over standard pressureless mat- 
ter |52| . The acceleration parameter becomes negative, 
according to current observations, when z < Zacc- More- 
over, the transition redshift, given by Eq. ([22|) . is in an 
excellent agreement with the value predicted by ACDM, 
i.e. Zacc, A — 0.75. The variation of q with respect to the 
cosmological redshift z, measures the rate of change of 
the acceleration, i.e. j{z). For our model, it is easy to 
get 

., ,_ Po + PiZ{z) + P^Z{zf + P:iZ{zf 

where for simplicity we defined the following constants 

Po = -Pa - 9a'(l - ^rr.f + 9a(l - fi™)'pA , 

/5i = (7 - 10nm)PA - 9a(l - Qm? , 

/52 = -1 - 2n^ , (24) 

P3 = ■ 

PA 

and the position 



Z(z) = a(l-an) (^ + 1) 



(25) 



In particular, j{z) is also known in literature as the jerk 
var am,eter and it is given by the definition j = -^ — 3q — 2 
[53| . We are interested in determining the jerk parameter 
at z = 0. It is expected that j > today in order to 




FIG. 1. Representation of supernova magnitude ij,{z) VS red- 
shift 2. Data and the best fit, both related to the correspond- 
ing uncertainties (red lines). The union 2.1 dataset has been 
plotted. 



guarantee that for z > Zacc the acceleration parameter is 
positive [5J, [55[ . The jerk parameter today reads 



pI - 9apAil - n.,n)^ + 9a2(i _ q,^^)3 
Pa 



(26) 



and the expected acceleration parameter today reads 
1 r 



qo 



Silr^ 



3T(l-r2„ 



1 



(27) 



By the definitions of go and jo, it is possible to find out 
viable priors for the free parameters, a, pA and flm- In 
particular, wc summarize the priors, that we impose in 
our computational analysis, in Tab. I, taking into ac- 
count the results obtained in Eqs. ^, ^ and (P7| . 



IV. COSMOLOGICAL CONSTRAINTS 

In this section, we constrain the free parameters of our 
model, through a numerical analysis based on current 
cosmological data. We rely on three statistical sets of 
parameters. In doing so, we define three different max- 
imum order of parameters, assuming a hierarchy among 
the three sets. The sets are summarized as follows 



A = 
B = 
C = 



Ho, ^m, 



PA 



PA 



PA 



(28) 



The hierarchy of Eqs. ([28|) predicts a broadening of 
the sampled distributions, adding the cosmological co- 
efficients to constrain. In the Gaussian regime, the error 
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FIG. 2. l-dimensional marginalized contour plots for Ho, -^ and flm, using A. We considered the union 2.1 compilation. Here, 
we report la, 2a and 3a of confidence levels. 



propagation becomes higher as hierarchy, between pa- 
rameters, increases. To alleviate the error propagation 
and possible systematics, we consider the numerical pri- 
ors reported in Tab. I. 



A. Priors on the cosmological parameters and the 
initial condition on Hq 

Viable cosmological priors are need in order to allevi- 
ate the so called degeneracy between cosmological dis- 
tances. In fact, the luminosity distance by itself is not 
enough to separately constrain all the cosmological den- 
sities. It follows that the EoS cannot be constrained with 
arbitrary accuracy with SNcIa data only. By fixing vi- 



able priors, one needs to reduce the total phase space. 
This permits us to complement distance measurements 
by different constraints. In other words, once the priors 
have been determined, it is possible to infer cosmologi- 
cal bounds from one class of measurements. Moreover, 
the simple choice of using different distances does not 
guarantee a priori that the physical region for constrain- 
ing the free parameters is actually reduced. For those 
reasons, we rely on three combined tests, performed by 
using SNela, BAO measurements and the CMB surveys. 
Our choice determines a combination of cosmological and 
geometrical procedures, which allows to circumscribe the 
phase space with higher precision. In addition, spatial 
geometry is also set to be geometrically flat, while the 
initial condition on Hq is determined as follows: 1) first, 
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FIG. 3. 1-dimensional marginalized contour plots for Hq, -^ and f2m, using set B. We considered the union 2.1 compilation. 
Here, we report la, 2cr and 3(t of confidence levels. 
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FIG. 4. We show the behavior of the EoS of our model, in terms of the redshift z (left figure). Moreover, we plot the Hubble 
rate (right figure), comparing it with the ACDM model. Respectively red line for our model and dashed line for ACDM. We 
used the indicative values Q.m — 0.290, -^ = 1.1 for our model, while for ACDM we used the values of WMAP 7-years. 



q(z) 




FIG. 5. Graphic of q{z) for our model (color line) and ACDM (dashed). We notice small differences in the redshift transitions 
from the acceleration to deceleration phases and good agreement with the numerical value at small redshift with respect to 
present cosmographic bounds. 



Flat priors 



0.5 < h < 0.8 

0.009 < fib/i^ < 0.08 

0.05 < fidmft^ < 0.26 

0.01 PA < Q < 100 PA 



Additional constraints 







0.33 



TABLE I. Priors imposed on the free parameters, involved 
in the the Bayesian analyses. We considered the additional 
geometrical assumption of a spatially flat universe. 



we consider Hq free to vary, 2) second, we constrained it 
by WMAP 7- years results, 3) finally, we impose its value, 
by assuming the Hubble space telescope (HST) measure. 
Additional cosmological and geometrical priors, adopted 
throughout our numerical analysis, are summarized in 
Tab. I, as already underlined. 



B. Cosmological datasets 

In our numerical analysis, we consider three different 
datasets. In particular, we take into account the union 
2.1 compilation, the BAO measure and the CMB surveys, 
with the constraint on Hq given by the measurement of 
the WMAP 7-years and HST respectively. In the union 
2.1 compilation of the supernova cosmology project [5a |. 
the covariance matrix has been evaluated with and with- 
out systematics. The compilation includes previous sur- 
veys, i.e. union 2 [STj and union 1 [5»|. The supernova 



measurements may be represented in the plane modulus- 
rcdshift, i.e. /i— z. They consist of 580 measurements of fi 
and z, spanning in the redshift range 0.015 < z < 1.414. 
The error over z is assumed to be negligibly small, while 
to each supernova is associated a corresponding error on 
/i. This turns out to be important since supernovas are 
standard indicators and represent the primary distance 
indicators. The relevance of supernova measurements is 
related to the fact that their rest frame wavelength region 
spans from 4000 to 6800 A for all the transient events 
which characterize supernovas. By assuming that the 
same rest frame wavelengths are measured at all z, one 
can compare the supernova brightness as independent as 
possible from a supernova model. Thence, to fix cosmo- 
logical constraints, we make use of a Bayesian method in 
which the best fits of parameters are inferred by maxi- 
mizing the following likelihood function 



C ex exp(-xV2) , 



(29) 



where x^ is the [pseudo)chi- squared function. To ob- 
tain the corresponding posterior distributions, we con- 
sider uniform priors in the interval ranges of Tab. I. It is 
easy to show that the luminosity distance is 



dL{z) = {l + z) 



m)' 



(30) 



where £ = j^, and by defining the distance modulus /i 
for each supernova 



^j = 25 + 5 log 



10 



Mpc ' 



(31) 



together with the corresponding l-cr^ error, we can 
rewrite the x^ parameter of Eq. (|29p as 



XSN 



E^ 



thcor 



bs\2 



t^i) 



(32) 



TABLE II. Best fits of the free parameters of our model, tested by SNela. The quoted errors show tlie 68.3%, 95.4% and 
99.7% confidence level uncertainties. 



Parameter 


Set A 
xLn = 0.9727 


Sets 
xLn = 0.9727 


Sets 
xLn = 0.9756 


SetC 
xLn = 0.9728 


SetC 
Xlin = 0.9837 


■Ho 


70 9 + 1.3 + 2.0 + 2.5 


Ho = 70.2 


Ho = 74.2 


-Ho = 70.2 


-Ho = 74.2 


fim 


n 97S + 0. 052+0. 083+0. 10 
"•^'"-0.048-0.072-0.088 


n 978+0. 052 + 0. 0S3 + 0. 10 
w.z, 1 u_o. 048-0. 072-0. OSS 


n '^n'i + 0. 038+0. 060 + 0. 075 
""^""^-0.036-0.055-0.068 


Qm = 0.274 


Q.m = 0.274 


a 

PA 


1 49 + 0.27+0.42 + 0.52 
-0.26-0.39-0.48 


1 4.9 + 0.27+0.41 + 0.52 
-0.26-0. 39-0. 48 


9 44+0.32+0.51 + 0.58 
-0.31-0.47-0.64 


1 40 + 0-26+0.41 + 0.51 
"- 0.25 -0.39 -0.47 


9 99 + 0-26+0-44+0.56 
-0.26-0.42-0.52 



Notes. "Ho is given in Kni/s/Mpc. 



TABLE III. Best fits of the free parameters with BAO. Note that with BAO it is not possible to estimate Ho- The quoted 
errors show the 68.3%, 95.4% and 99.7% confidence level uncertainties. 



Parameter 


Set y^ 


SetC 


fi™ 


n 98^^ + 0.030+0.043+0.053 
"• ^^-0.027-0.038-0.046 


n„ = 0.274 


P^ 


1 QQ + 0. 36+0. 59+0. 74 
■-^^-0.41-0.50-0.59 


1 48+0.42+0.61+0.76 
" -0.37-0. 51-0. 62 



A simple test with SNcIa minimizes the term Xsn^ which 
corresponds to maximize Eq. ((29|) . 

The large scale galaxy clustering observations provide 
the signatures of the BAO [5^. The theoretical back- 
ground employs that the universe consisted of a hot 
plasma of photons, electrons, protons, baryons and other 
light nuclei, at a certain epoch of its evolution. As a con- 
sequence, the Thompson scattering between photons and 
electrons leads to oscillations in the hot plasma. As the 
universe becomes neutral, it is possible to consider the 
initial perturbation patterns which are imprinted on the 
matter distribution. Hence, by observing the spectrum 
of galaxy correlations today, it is possible to focus on 
observations of BAO. The corresponding measurement is 
represented by a sound horizon length, whose physical 
meaning deals with the distance traveled by an acous- 
tic wave by the time of plasma recombination. The BAO 
peak is considered a standard cosmological ruler, because 
such a peak is independent of the choice of a particular 
cosmological model. We use the peak measurement of 
luminous red galaxies, denoted by Abao- It reads 

,1 r 1 
Abao ={^b + ^dti 



1 



^izBAO 

AO -1 



d^ 



(33) 



.zbao Jo ^iO 

with fib and fldm respectively the baryonic and dark mat- 
ter densities and zbao = 0.35. In addition, the observed 
Abao is estimated to be 

ABAo,ots = 0.469 ij^] = 0.477 , (34) 



^0.987 

with an error (ta ~ 0.017. In the case of the BAO mea- 
surement, we minimize the chi square 

2 



Xbao — 



A — Aq, 
(TA 



bs 



(35) 



The third test is represented by CMB observations. This 
kind of measurements has recently reached much interest, 
as a nearly isotropic background was discovered in 1965. 
Measures of the WMAP satellite found new data, able 
to alleviate the cosmological degeneracy between mod- 
els. The underlying philosophy relies on CMB radiation, 
which can be directly detected by keeping in mind that 
it is influenced by two cosmological epoches, i.e. the last 
scattering era and present time [60, |6l[ . For the CMB 
test, we define the so-called CMB shift parameter Rcmb- 
Its standard definition reads 



Rcmb = {VLb + ^dm) ' 



£(0 



(36) 



with Zrec ~ 1000. Equation (|36|) presents several the- 
oretical shortcomings [62|. It is interesting to replace 

such a relation with Rcmb = 2+, where li is the posi- 

'i 
tion of the first peak on the CMB TT power spectrum of 

the model under consideration. Moreover, li is the first 

peak in a flat homogeneous and isotropic universe with 

flf) = I — ^drn- Hereafter, li is written as 

h = DA{Zrec)s{Zrec)~^ , (37) 

where DA{zrec) is called co- moving angular distance, i.e. 



DAiZrec)^ / (l+O^C, 

Jo 



(38) 



with s{zrec) representing the sound horizon at recombi- 
nation 



1 f°° 



(39) 



We defined o(^) as the sound speed of the photon-to- 
baryon fluid. Its expression reads D = f 3 + 4y- 



TABLE IV. Best fits of tlie free parameters witli CMB. Note that witli CMB it is not possible to estimate Tio- The quoted 
errors show the 68.3%, 95.4% and 99.7% confidence level uncertainties. 



Parameter 


Set A 


SetC 


fim 


D 9R1 +0.034+0.050+0.062 
^^^ -0.030-0.041-0.050 


fi™ = 0.274 


a 


1 1 7+0.30+0.45 + 0.57 
' -0.23-0.31-0.37 


1 94+0.31 + 0.47+0.60 
-0.24-0.33-0.39 



TABLE V. Best fits of the parameters for the three considered models with SNela, BAO and CMB. The quoted errors show 
the 68.3%, 95.4% and 99.7% confidence level uncertainties. 



Parameter 


Set A 
xLn = 0.9727 


Set B 
xLn = 0.9727 


Sets 
xLn = 0.9871 


SetC 
xLn = 0.9729 


SetC 
Xlin = 0.9885 


no 


70 2 + 1-3 + 2.0 + 2.5 
' ^ -1.3-2.0-2.4 


Ho = 70.2 


Ho = 74.2 


-Ho = 70.3 


-Ho = 74.2 


n™ 


978 + 0-049+0.078 + 0.099 
^ ' ^-0.045-0.068-0.083 


n 97R+0. 049 + 0. 07S + 0. 99 
^' ' ^-0.045-0.068-0.083 


fl 9SC; + 0. 032+0. 056+0. 072 
^' ^^-0.031-0.051-0.064 


Qm = 0.274 


n„ = 0.274 


PA 


1 4.9 + 0.26 + 0.41 + 0.51 
-0.25-0.39-0.48 


1 4.9 + 0.26 + 0.41 + 0.51 
-0.25-0.39-0.48 


9 9S + 0. 25+0. 44 + 0. 57 
"- 0.25 -0.42 -0.53 


1 4^0 + 0.26+0-40+0.51 
^ ^-0.25-0.38-0.47 


O 01+0.24+0-43+0.55 
-0.23-0.40-0.51 



Notes. "Ho is given in Kni/s/Mpc. 



The importance of using CMB leads to its complemen- 
tary with respect to SNela and BAO measurements. This 
comes from the fact that SNela and BAO are confined 
in low rcdshift regimes, constrained to z < 2, while for 
CMB, the recombination redshift is higher of three orders 
of magnitude |63|. According to the CMB measurement, 
we minimize the following 



2 / /v I'Cobs 

XCMB 



(40) 



Since the three different sets of observations are not cor- 
related between them, the total x^ reads 



X 



i=l 



2 2 2 

XSN + XbAO + XcMB 



(41) 



From one hand, the SNcIa test needs to fix _ffo, as ini- 
tial condition, in addition to the free parameters of our 
model, i.e. a and p\. It is prominent to assume gaus- 
sian priors on the Hubble constant, in order to fulfill 
the initial condition on Hq. In doing so, we consider 
the WMAP 7-years measurements, which suggests "Hp = 
70.2±3.6fcm/s/Mpcandifo = 74.2±3.6 fcm/s/Mpc [ll], 
as measured by the HST. Following these constraints, and 
those of Tab. I, wc perform three separate tests with 
SNela, BAO and finally CMB, by combining them to- 
gether. On the other hand, the cosmological tests BAO 
and CMB do not directly depend on Hq. It follows that 
the unique bound on Hq should be imposed when we 
perform the supernova test only. 



V. MODEL COMPARISON THROUGH AIC 
AND BIC SELECTION CRITERIA 

As already stressed in Sec. Ill, the degeneracy prob- 
lem, between cosmological models, plagues standard 



techniques of numerical analyses. The task of discrim- 
inating which models better fit current data is actually a 
thorny issue of statistics. The main disadvantage lies on 
assuming, in a numerical computation, a particular form 
oi H{z) for the cosmological model. In other words, each 
cosmological test assumes a priori that the cosmologi- 
cal model under exam is statistically favored. This is 
a consequence of the y^ analysis, which is able to con- 
strain the free parameters of a given model, although it 
does not provide any information about the validity of 
the model itself. To alleviate this problem, one could 
use model independent procedures of statistical analysis. 
Among various possibilities, there exist in the literature 
statistical methods, able to understand which model is 
really favorite, than others, once the same data survey 
is used. It is easy to show that comparing different chi 
squares among alternative models, with the same num- 
bers of free parameters, it is possible to check which 
model is statistically favored than others. In lieu of lim- 
iting our attention to simply fit our model, as in Sec. IV, 
wc compare our approach with two relevant paradigms, 
i.e. a variable quintessence model: the Chevallier, Po- 
larsky, Linder (CPL) parametrization [7^, UM, ) and the 
standard ACDM model. Naively, we can expect that a 
combination of lowest chi squares and fewest numbers of 
parameters provides the model which better reproduces 
cosmological data. We summarize below the ACDM and 
CPL Hubble rates respectively, i.e. 



£i 



Q.„,{i + zf + 1 ~ n^ 



^r.Al + zf + {l-n^)i{z) 



(42) 



with f(z) = (1 + z)3{i+^o+^a) exp (-iuja-p^\ . In partic- 
ular, oj — —1 and w — uq + ujai jr-^ ) , respectively for 
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TABLE VI. BIG and AIC analysis, performed by assuming model A with SNela data. 
Model Num. of Par. (k) Parameters Xmm ABIC AAIC 



ACDM 

Ising Fluid 

GPL 



1 Qm 0.9727 

2 !r2„, ^ 0.9727 6.363 2 

3 fi™, ivo, u}a 0.9879 12.741 4.015 



ACDM and CPL models. In what follows, wc describe 
the statistical methods that we are going to use, in order 
to make a comparison between our model and the above 
cited ACDM and CPL paradigms. 



A. Selection criteria 

In this subsection, we propose two statistical methods 
which follow the guidelines given in Sec. V. They are 
the so-called AIC [H^ and BIC [^ selection crite- 
ria. The first one has reached a widely accented con- 
sensus, becoming a common diagnostic tool l70l473l . It 
has been largely used for regression models [7J|, since 
its first applications |65l - [68| . The basic demands of AIC 
and BIC consist in postulating two statistical distribu- 
tions, i.e. /(x) and g{x\6). The first distribution f{x) 
is imposed to be the exact reconstruction of a particu- 
lar model, while the second distribution, i.e. g{x\9), is 
thought to approximate /(a;), by bounding A'^ parame- 
ters, that are included within the vector 9. Once f{x) 
and g{x\9) are defined by numerical procedures, the set 
of parameters 9 is estimated, minimizing the departures 
between f{x) and g(x,9min)- 

Obviously, the function f{x) is not known a priori. 
Hence, both the AIC and BIC criteria are meaning- 
less if evaluated for a single model. In other words, 
once the minima of AIC and BIC are determined, i.e. 
AlCmin and BICmin respectively, we should investigate 
the differences AA/C = AIC - AlCmin and ABIC = 
BIC — BICmin ■ A general form for the AIC function is 

AIC = -21nC,na. + 2k, (43) 

and by following [63, we define the BIC function as 



BIC 



-21n/:^ 



klnN , 



(44) 



where for both the methods AIC and BIC, Cmax is the 
maximum likelihood function, corresponding to the min- 
imum of x^. Moreover, k is the number of parameters 
to be estimated, and N is the number of data points 
used to perform the cosmological fits. To Gaussian 
distributed errors corresponds Xmin ^ ~'^^"^^max and 
we can simphfy AAIC = Ax^„j -I- 2Afc and ABIC = 
AxLn + AfchrAT. 

In Tab. Ill, wc report the numerical results of using the 
AIC and BIC criteria. In particular, we find for Ei and 
£2 respectively: f7,„ = 0.2791°;°^^, Ho = mMQt\fm, 



0.9727, a„, = 0.282+^:^?^^ Ho = 70.120+1:^1^ 



Wo 



-0.95 



+0.49 , , 

^0.46' ^a 



1.52+i:^^ and xsn, 



0.9879. 



We find that our model fairly good adapts to cosmo- 
logical data, being disfavored with respect to the stan- 
dard ACDM picture, and behaving better than CPL. 
The statistical success of ACDM is clearly due to the 
smallest number of parameters involved into calcula- 
tions. However, our paradigm shows small departures 
from ACDM, behaving much better than an evolving 
quintessence model. This can be interpreted as a pos- 
sible indication that our network interacting Ising fluid 
may be seen as a relevant alternative to ACDM. 



VI. FINAL REMARKS 

In this work, we investigated the possibility to model 
DE through an Ising fluid on a lattice with network in- 
teractions. A negative pressure, associated to DE and 
compatible with present observations, emerged by con- 
sidering the barotropic EoS of our model in the equilib- 
rium configuration. We inferred theoretical constraints in 
a fairly good agreement with current observations. We 
demonstrated that, at low redshift, it is possible to obtain 
a viable Hubble rate which well mimics the DE effects, re- 
ducing to ACDM at the zero order expansion, in terms of 
the DE density. The corresponding acceleration parame- 
ter and its variation, namely the jerk parameter, provided 
interesting results which properly fitted with modern ob- 
servations. In particular, the acceleration starts at a red- 
shift which is excellently close to the one predicted by 
ACDM. In addition, the present value of j(z) confirmed 
that the acceleration parameter has changed its sign, as 
the universe expands. From such considerations, it fol- 
lowed that one of the main advantages of our model re- 
lied on interpreting the DE nature as an emergent Ising 
fluid with network interactions. In addition, we fixed 
constraints on the free parameters of our model, by em- 
ploying three cosmological datasets, i.e. the SNela, BAO 
and CMB surveys. In particular, we first used SNela, 
BAO and CMB separately and then we combined SNela 
with BAO and CMB, constraining the free parameters of 
our model in tighter intervals, though viable geometri- 
cal and cosmological priors. In doing so, we chose for the 
computational analysis three cosmological sets of observ- 
ables, ordered in a hierarchial way, evaluating the corre- 
sponding errors up to 3a of confidence level. The results 
confirmed that our model could be viewed as a viable al- 
ternative to ACDM, in order to describe the DE effects at 
present time. To this end, we showed that, through the 
use of the AIC and BIC selection criteria, our model pro- 
vided small departures than ACDM, behaving smoother 
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than the so called CPL parametrization. Future efforts 
could be devoted to apply the network interacting Ising 



model to other stages of the universe evolution, wonder- 
ing whether the model could be modeled for different 
epochs of the universe evolution. 
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